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Abstract 



We present a quantitative study of the ability of models with different levels of 
hierarchy to reproduce the solar neutrino solutions, in particular the LA solu- 
Q_i! tion. As a flexible testing ground we consider models based on SU(5)xU(l)p. 

In this context, we have made statistical simulations of models with different 



patterns from anarchy to various types of hierachy: normal hierarchical mod- 
els with and without automatic suppression of the 23 (sub) determinant and 
inverse hierarchy models. We find that, not only for the LOW or VO solu- 
tions, but even in the LA case, the hierarchical models have a significantly 
better success rate than those based on anarchy. The normal hierachy and 
the inverse hierarchy models have comparable performances in models with 
see-saw dominance, while the inverse hierarchy models are particularly good 
in the no see-saw versions. As a possible distinction between these categories 
of models, the inverse hierarchy models favour a maximal solar mixing angle 
and their rate of success drops dramatically as the mixing angle decreases, 
while normal hierarchy models are far more stable in this respect. 
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1 Introduction 



At present there are many possible models of v masses and mixing [1 . This variety is 
mostly due to the considerable experimental ambiguities that remain. In particular 
different solutions for solar neutrino oscillations are still possible. Although the 
Large Angle (LA) solution emerges as the most likely from present data, the other 
solutions LOW and Vacuum Oscillations (VO) are still not excluded (the small angle 
solution is very disfavoured by now and we will disregard it in most of the following 
discussion). Indeed no solution is actually leading to an imposingly good fit and, 
for example, the discrimination between LA and LOW is only based on a few hints 
which are far from compelling J2] . Hopefully in a few months, when the first results 
from the KamLAND experiment [3| will be known, one will have decisive evidence 
on this matter. Here we tentatively disregard the possibility of a third neutrino 
oscillation frequency as indicated by the LSND experiment [I] but not confirmed 
by KARMEN j3] and to be finally checked by the MiniBOONE experiment jB] now 
about to start. 

For model building there is an important quantitative difference between the LA 
solution on the one side and the LOW or VO solutions on the other side. While for 
all these solutions the relevant mixing angle # 12 is large, the value of the squared 
mass difference AmJ 2 = m\ — m\ (with, by definition, m 2 > mf) is very different for 
LA, LOW and VO: ~ KT 4 eV 2 , ~ KT 7 eV 2 and ~ KT 10 eV 2 , respectively. Thus the 
gap with respect to the atmospheric neutrino oscillation frequency Am^ = m\ — m\, 
which is given by |Am| 3 | ~ 3 ■ 10~ 3 eV 2 , is moderate for LA and very pronounced 
for the other two solutions. 

For the LOW and VO solutions the large frequency difference with respect to that 
of atmospheric neutrinos points to a hierarchical spectrum for the three light neu- 
trinos. Possible hierarchical patterns are the normal hierarchy case ^> ni2 ^ m\ 
or the inverted hierarchy alternative m 2 ^ mi > m 3 (in this case Am| 3 is negative 
in our definition). Then a main problem is to explain the presence of large mixing 
angles between largely splitted mass states (in particular the atmospheric neutrino 
oscillation mixing angle #23 is experimentally close to maximal). In hierarchical mod- 
els the consistency of these usually opposed constraints is obtained by mechanisms 
that guarantee a vanishing or a strongly suppressed 23 sub-determinant. In the 
see-saw mechanism for neutrino masses this suppression can be naturally obtained, 
for example, in the so-called lopsided models [7j and/or by the dominance jH] of 
one eigenvalue in M~ l , M being the right-handed (RH) Majorana matrix. Models 
of this type have been studied and provide, as also quantitatively confirmed by our 
present analysis, an essentially unique framework for a successful description of both 
atmospheric and solar neutrino oscillations when the LOW or the VO solutions are 
adopted. 



1 



In the case of the LA solution the ratio of the solar and atmospheric Am 2 ranges 
is typically given by 

r = Am? 2 / 1 Am 2 23 \ ~ 1/20- 1/100. (1) 

For LA one can reproduce the data either in a nearly degenerate or in a hierarchical 
model. In a degenerate model, due to laboratory and cosmological bounds, the 
common value of m 2 ~ m 2 cannot exceed a few eV 2 . But the actual value is 
probably well below this level because of the constraint imposed by neutrinoless 
double beta decay (0^/3/3) [5] that would otherwise require a strong cancellation, 
only possible for nearly maximal solar oscillation mixing jTU]. This fact, together 
with the general difficulty, in the absence of a specific mechanism, of obtaining 
too small values of Am 2 2 /m 2 , suggests that a moderate degeneracy is more likely. 
Typically we could have all mf ~ (few 10~ 3 — 10~ 2 ) eV 2 with one comparatively 
not-so-small splitting Am 2 2 ~ 10~ 4 eV 2 . Or, as a different example, we can have a 
(normal) hierarchical model with m 2 ~ 3 ■ 10~ 3 eV 2 , m 2 . ^ 10~ 4 eV 2 and m\ ~ or 
an (inverse) hierarchical model with m 2 2 ~ 3 ■ 10~ 3 eV 2 , and m 2 ~ 0. Actually a 
sufficient hierarchy (a factor of 5 in mass) can arise from no significant underlying 
structure at all. In particular, the see-saw mechanism, being quadratic in the Dirac 
neutrino masses, tends to enhance small fluctuations in the Dirac eigenvalue ratios. 
This is the point of view of anarchical models [TT], where no structure is assumed 
to exist in the neutrino mass matrix and the smallness of r is interpreted as a 
fluctuation. But one additional feature of the data plays an important role in this 
context and presents a clear difficulty for anarchical models. This is the experimental 
result that the third mixing angle #13 is small, |{7 e3 | = | sin #13! ^ 0.2 [12 . So, for 
neutrinos two mixing angles are large and one is small. Instead in anarchical models 
all angles should apriori be comparable and not particularly small. Therefore this is 
a difficulty for anarchy and, for the survival of these models, it is necessary that # 13 
is found very close to the present upper bound. Instead in hierarchical models the 
smallness of #13 can be obtained as a reflection of the underlying structure in that 
some small parameter is present from the beginning in these models. 

In this note we make a quantitative study of the ability of different models to 
reproduce the solar neutrino solutions. As a flexible testing ground we consider 
models based on SU(5)xU(1)f- The SU(5) generators act "vertically" inside one 
generation, while the U(1)f charges are different "horizontally" from one generation 
to the other. If, for a given interaction vertex, the U(1)f charges do not add to zero, 
the vertex is forbidden in the symmetric limit. But the symmetry is spontaneously 
broken by the VEV v / of a number of "flavon" fields with non vanishing charge. Then 
a forbidden coupling is rescued but is suppressed by powers of the small parameters 
Xf = Vf/A with the exponent larger for larger charge mismatch We expect 

Vf ft Mqut an d, for the cut-off A of the theory, A ^ M Pi . In these models [HJ [TS] 
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the known generations of quarks and leptons are contained in triplets and , 
(i = 1,2,3) corresponding to the 3 generations, transforming as 10 and 5 of SU(5), 
respectively. Three more SU(5) singlets describe the RH neutrinos. In SUSY 
models we have two Higgs multiplets, which transform as 5 and 5 in the minimal 
model. All mass matrix elements are of the form of a power of a suppression factor 
times a number of order unity, so that only their order of suppression is defined. We 
restrict for simplicity to integral charges: this is practically a forced choice for the 
LA case where the hierarchy parameter must be relatively large (so that v^A ~ 0(1)), 
while for the LOW and VO cases, where the hierarchy parameter is small, it is only 
motivated by the fact that enough flexibility is obtained for the present indicative 
purposes. There are many variants of these models PQ: fermion charges can all be 
non negative with only negatively charged flavons, or there can be fermion charges 
of different signs with either flavons of both charges or only flavons of one charge. 
The Higgs charges can be equal, in particular both vanishing or can be different. 
We will make use of this flexibility in order to study the relative merits of anarchy 
versus various degrees and different patterns of hierarchy. 

In this context we have studied in detail different classes of models: normal hier- 
archical models with and without automatic suppression of the 23 (sub) determinant, 
inverse hierarchy models and anarchical models. The normal hierarchical models 
without automatic suppression of the 23 determinant are clearly intermediate: in a 
sense in those cases anarchy is limited to the 23 sector. We denote them as partially 
hierarchical or semi-anarchical in the following. We also compare, when applicable, 
models with light neutrino masses dominated by the see-saw mechanism or by non 
renormalizable dim-5 operators. We construct our models by assigning suitable sets 
of charges for \E^°, and ^j. In all input mass matrices the coefficients multiplying 
the power of the hierarchy parameter are generated at random as real and complex 
numbers in a given range of values fTJ • We compare the case of real or complex 
parameters and we also discuss the delicate questions of the probability distribution 
for the coefficients and the stability of the results. We assign a merit factor to each 
model given by the percentage of success over a large sample of trials. For each 
model the value of the hierarchy parameter is adjusted by a coarse fitting procedure 
to get the best rate of success. 

Our results can be summarized as follows. As expected, for the LOW and VO 
cases only hierarchical model provide a viable approach: in comparison the rate 
of success of anarchical and seminarchical models is negligible. But also for the 
LA solution we still find that hierarchical models are sizeably better in general. 
The most efficient ones are inverse hierarchy models with no see-saw dominance, 
which are more than 10 times better with respect to anarchy with see-saw (anarchy 
prefers the see-saw case by about a factor of 2). Among the see-saw dominance 
versions the most performant models remain the hierarchical ones (by a factor of 
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about 4 with respect to anarchy with see-saw) with not much difference between 
inverse or normal hierarchy. Semi- Anarchical models are down by a factor of about 
2 with respect to hierarchical models among the see-saw versions (but this value is 
less stable with respect to changes of the extraction procedure and, for example, it 
tends to be washed out going from complex to real coefficients). In all models the 
tan 2 #23 distribution is in agreement with large mixing but it is not sharply peaked 
around 1 as for maximal mixing. Near maximal mixing is instead a prediction for 
solar neutrinos in inverse hierarchical models, so that their advantage with respect 
to other models would be rapidly destroyed if the data will eventually move in the 
direction away from maximality. 



2 Framework 

We consider a class of models with an abelian flavour symmetry compatible with 
SU(5) grand unification. Here we will not address the well-known problems of grand 
unified theories, such as the doublet-triplet splitting, the proton lifetime, the gauge 
coupling unification beyond leading order and the wrong mass relations for charged 
fermions of the first two generations. We adopt the SU(5)xU(1)f framework simply 
as a convenient testing ground for different neutrino mass scenarios. In all the 
models that we study the large atmospheric mixing angle is described by assigning 
equal flavour charge to muon and tau neutrinos and their weak SU(2) partners (all 
belonging to the 5 = (l,d c ) representation of SU(5)). Instead, the solar neutrino 
oscillations can be obtained with different, inequivalent charge assignments and both 
the LOW 1 and the LA solution can be reproduced. 

A first class of models is characterized by all matter fields having flavour charges 
of one sign, for example all non negative. An important property of models in this 
class is that the light neutrino mass matrix L T m v L is independent from the charges 
of both the 10 = (q, u c , e c ) and 1 = u c representations , even in the see-saw case when 
m u = nr^M^mn- For m„ entries the powers of the symmetry breaking parameter 
A are dictated by the charges F of the 5. Since in this case what really matters are 
charge differences, rather than absolute values, the equal charges for the second and 
third generations can be put to zero, without loosing generality: 

F(5) = (6,0,0) b>0 . (2) 

If b also vanishes, then the light neutrino mass matrix will be structure-less and 
we will call anarchical (A) this sub-class of models. In a large sample of anarchical 

1 The LOW and VO solutions can be fitted almost equally well by suitable models. Thus here 
we will focus mainly on the LOW case. 
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models, generated with random coefficients, the resulting neutrino mass spectrum 
can exhibit either normal or inverse hierarchy. Anarchical models clearly prefer the 
LA solution with a moderate separation between atmospheric and solar frequencies. 
They tend to favour large, not necessarily maximal, mixing angles, including U e 3, 
which represents a problem. Therefore, in anarchical models, U e s is expected to be 
close to the present experimental bound. 

If b is positive, then the light neutrino mass matrix will be structure-less only in 
the (2,3) sub-sector and we will call semi-anarchical (SA) the corresponding models. 
In this case, the neutrino mass spectrum has normal hierarchy. However, unless the 
(2,3) sub-determinant is accidentally suppressed, atmospheric and solar oscillation 
frequencies are expected to be of the same order and, in addition, the preferred 
solar mixing angle is small. Nevertheless, such a suppression can occur in a fraction 
of semi-anarchical models generated with random, order one coefficients. The real 
advantage over the fully anarchical scheme is represented by the suppression in U^. 

In a second class of models matter fields have both positive and negative flavour 
charges. In these models, the light neutrino mass matrix will in general depend 
also on the charges of 10 and 1. A first sub-case arises when only the RH neutrino 
fields have charges of both signs. It has been shown that it is possible to exploit 
this feature to reproduce a neutrino mass spectrum with normal hierarchy and a 
natural gap between atmospheric and solar frequencies. Via the see-saw mechanism 
the (2,3) sub-determinant vanishes in the flavour symmetric limit. At the same time 
a large solar mixing angle can be obtained. Clearly this is particularly relevant for 
the LOW (or VO) solution. It is less clear to which extent the condition of vanishing 
determinant is needed for the LA solution and one of the purposes of the present 
paper is precisely to compare these models, which we will call hierarchical (H) with 
the anarchical and semi-anarchical models, that do not reproduce such a condition. 

Finally, we can have fields with charges of both signs in both the 1 and the 5. In 
this context it is possible to reproduce an inverse hierarchical spectrum, with a large 
(actually, almost maximal) solar mixing angle and a large gap between atmospheric 
and solar frequencies [T2j. Also this sub-class of models, which we call inversely 
hierarchical (IH), are appropriate to describe both LOW and LA solutions. 

The hierarchical and the inversely hierarchical models may come into several 
varieties depending on the number and the charge of the flavour symmetry breaking 
(FSB) parameters. Here we will consider both the case (I) of a single, negatively 
charged flavon, with symmetry breaking parameter A or that of two (II) oppositely 
charged flavons with symmetry breaking parameters A and A'. In case I, it is impossi- 
ble to compensate negative F charges in the Yukawa couplings and the corresponding 
entries in the neutrino mass matrices vanish. Eventually these zeroes are filled by 
small contributions, arising, for instance, from the diagonalization of the charged 
lepton sector or from the transformations needed to make the kinetic terms canon- 
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Table 1: Models and their flavour charges. 



Model 


#10 


#5 


#i 


(H u , H d ) 


Anarchical (A) 


(3,2,0) 


(0,0,0) 


(0,0,0) 


(0,0) 


Semi-Anarchical (SA) 


(2,1,0) 


(1,0,0) 


(2,1,0) 


(0,0) 


Hierarchical (H) 


(3,2,0) 


(2,0,0) 


(i,-i,o) 


(0,0) 


Inversely Hierarchical (LA) 


(3,2,0) 


(1,-1,-1) 


(-i,+i,o) 


(o,+i) 


Inversely Hierarchical (LOW) 


(2,1,0) 


(2,-2,-2) 


(-2,+2,0) 


(0,+2) 



ical. In our analysis we will always include effects coming from the charged lepton 
sector, whereas we will neglect those coming from non-canonical kinetic terms. 

Another important ingredient in our analysis is represented by the see-saw mech- 
anism ^Hl- Hierarchical models and semi-anarchical models have similar charges in 
the (10, 5) sectors and, in the absence of the see-saw mechanism, they would give rise 
to similar results. Even when the results are expected to be independent from the 
charges of the RH neutrinos, as it is the case for the anarchical and semi-anarchical 
models, the see-saw mechanism can induce some sizeable effect in a statistical analy- 
sis. For this reason, for each type of model, but the hierarchical ones (the mechanism 
for the 23 sub-determinant suppression is in fact based on the see-saw mechanism), 
we will separately study the case where RH neutrinos are present and the case where 
they are absent. When RH neutrinos are present, there are two independent con- 
tributions to the light neutrino mass matrix. One of them comes via the see-saw 
mechanism from the exchange of the heavy RH modes. The other one is provided by 
L-violating dimension five operators arising from physics beyond the cut-off. These 
contributions have the same transformation properties under the flavour group and, 
in general, add coherently. In our analysis we will analyze the case where the see- 
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saw contribution is the dominant one (SS). The absence of RH neutrinos describes 
the opposite case, when the mass matrix is saturated by the non-renormalizable 
contribution (NOSS). 

For each type of model we have selected what we consider to be a typical rep- 
resentative 2 and we have collected in table 1 the corresponding charges. In the 
next section we will compare the performances of the following models: Ass, Anoss, 
SAss, SA N oss, H(ss,i), H(ss,ii), IH(ss,i), IH(ss,ii), IH(noss,i) and IH( NO ss,ii)- 

Table 2: Order of magnitude predictions for oscillation parameters, from neutrino 
mass matrices in eq. (j3J) and (jU) |lj d 2 3 denotes the sub-determinant in the 23 
sector and we show the effect of its accidental suppression for the semi-anarchical 
model. In the estimates we have chosen A = A'. Inverse hierarchy predicts an almost 
maximal (9 12 . 



Model 


parameters 


d 2 3 


Am 2 2 /|Am 2 3 | 


u e3 


tan 2 6 12 


tan 2 8 23 


A 


e = 1 


0(1) 


0(1) 


0(1) 


0(1) 


0(1) 


SA 


e = A 


0(1) 


0(4,) 


0(A) 


0(A 2 A4) 


0(1) 


Hn 


e = A 2 


0(A 2 ) 


0(A 4 ) 


0(A 2 ) 


0(1) 


0(1) 


Hj 


e = A 2 





0(A 6 ) 


0(A 2 ) 


0(1) 


0(1) 


IH (LA) 


e = r] = A 


0(A 4 ) 


0(A 2 ) 


0(A 2 ) 


1+0(A 2 ) 


0(1) 


IH (LOW) 


e = r] = A 2 


0(A 8 ) 


0(A 4 ) 


0(A 4 ) 


1+0(A 4 ) 


0(1) 



2 We made no real optimization effort to pick up the 'most' representative, but rather a model 
with a high success rate in its class. 



7 



Anarchical, semi-anarchical and hierarchical models give rise to a mass matrix 
for light neutrinos of the type 



V 



6 e\ 
1 1 
1 1 



(A, SA,H) 



(3) 



J 



where all the entries are specified up to order one coefficients and the overall mass 
scale has been conventionally set to one. For anarchical models, e = 1. Then 
all the entries are uncorrelated numbers of order one and no particular pattern 
becomes manifest. For the semi-anarchical model of table 1, e = A. There is a clear 
distinction between the first row and column and the 23 block of the mass matrix, 
which is structureless as in the anarchical models. In particular, barring accidental 
cancellations, the sub-determinant in the 23 sector is of order one. Finally, the 
hierarchical model defined by the choice of charges in table 1, has e = A 2 . At 
variance with the anarchical or semi-anarchical models, the determinant of the 23 
sector is suppressed by the see-saw mechanism and is of order AA'. 

The inversely hierarchical models are characterized by a neutrino mass matrix 
of the kind 



e 2 1 



7] 7] 



1 \ 

2 



(IH) 



(4) 



y 1 tf tf ) 



where e = A (A 2 ) and rj — \' (A /2 ) for the LA (LOW) solution. The ratio between 
the solar and atmospheric oscillation frequencies is not directly related to the sub- 
determinant of the block 23, in this case. The above mass matrices also receive an 
additional contribution from the diagonalization of the charged lepton sector, which, 
however, does not spoil the displayed structure. For completeness, we collect in table 
2 the gross features of the models under consideration. Notice that the hierarchical 
models predict a ratio Am 2 2 /| Amiy of order A 4 or A 6 . In these cases it is possible 
to fit both the LOW and the LA solutions with an expansion parameter A A' that, 
within a factor of two, matches the Cabibbo angle. On the contrary, the inversely 
hierarchical model that reproduces the LA solution favours Am 2 2 /|Am 2 3 | « 0(A 2 ). 
If we adopt this model to fit also the LOW solution, the corresponding values of 
A ~ A' do not provide a decent description of the remaining fermion masses. For 
this reason, when analyzing the LOW solution, we have considered a separate set of 
charges for the inverted hierarchy. 
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3 Method and results 



Abelian flavour symmetries predict each entry of fermion mass matrices up to an 
unknown dimensionless coefficient. These coefficients, that are free-parameters of 
the theory, are expected to have absolute values of order one. Aside from generalized 
kinetic terms, which we do not consider, the relevant mass matrices are specified 
by iV = 24 order one parameters: 9 from the charged lepton sector, 9 from Dirac 
neutrino mass matrix entries and 6 from Majorana RH neutrino mass matrices. 
When RH neutrinos are absent (that is, for the NOSS cases) the LH neutrino mass 
matrix contains 6 parameters and the relevant set reduces to N = 15 parameters. 
We have analyzed the case of real or complex coefficients V = (Vi, ...Vn), with 
absolute values generated as random numbers in an interval X and random phases 
taken in [0,27r]. To study the dependence of our results on X, we have considered 
several possibilities: [0.5,2] (default), [0.8,1.2], [0.95,1.05] and [0,1]. In the case 
of real coefficients, which is studied for comparison, we allow both signs for the 
coefficients. For each model, only a portion V$ °f the volume V = (2ttX) n of the 
parameter space gives rise to predictions in agreement with the experimental data 
within the existing uncertainties. We may interpret Vs/V as the success rate of the 
model in describing neutrino data. Clearly this portion shrinks to zero for infinitely 
good measurements. Therefore we are not interested in its absolute size, but rather 
in the relative sizes of Vs/V in different models. 

We evaluate the success rate of each model by considering, through a random 
generation, a large number of 'points' V and by checking whether the corresponding 
predictions do or do not fall in the experimentally allowed regions ^E|- To this pur- 
pose we perform a test based on four observable quantities: Ox — r = Am^/I Amjy, 
2 = tan 2 6> 12 , 03 = \U e3 \ = | sin#i 3 | and 4 = tan 2 6> 23 . We take [T§] : 



4- 10" 6 < r < 4 - 10~ 5 
\U e3 \ < 0.2 

0.52 < tan 2 fl 12 < 1.17 
0.33 < tan 2 6 23 < 3.3 



(LOW) (5) 



(LA) (6) 



0.01 < r < 0.2 

\U e3 \< 0.2 

0.24 < tan 2 fl 12 < 0.89 
0.33 < tan 2 ^3 < 3.3 

The boundaries of these windows are close to the 3a limits on the corresponding 
observable quantity. The test is successfully passed if Oi(V) are in the above win- 
dows, for all i = 1, ...,4. We will study the dependence on the choice made in eqs. 
(|5I6J) . by analyzing the distributions of the points generated for each observable. We 
then estimate Vs/V from the ratio between the number of successful trials over the 
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total number of attempts. Notice that we do not extend our test to the remaining 
fermion masses and mixing angles, though the models considered can also repro- 
duce, at the level of order of magnitudes, charged lepton masses, quark masses and 
mixing angles. With the interval X fixed at its reference value, [0.5, 2], the different 
models are compared at their best performance, after optimizing for each model the 
value of the symmetry breaking parameters A and A'. 



3.1 LOW solution 

In fig. we compare the success rates for the LOW solution. In this figure the 
anarchical or semi-anarchical models do not appear simply because their rates of 
success are negligible on the scale of the figure. 
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IH2 
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IH2 



IH1 




Figure 1: Relative success rates for the LOW solution. The sum of the rates has 
been normalized to 100. The results correspond to the default choice X = [0.5,2], 
and to the following values of A = A': 0.1, 0.15, 0.03, 0.04, 0.05, 0.06 for the models 
H(ss,ii), H (ss ,i), IH( SS ,ii), IH ( ss,i), IH (NO ss,ii) and IH (NOS s,i), respectively. The error 
bars represent the linear sum of the systematic error due to the choice of X and the 
statistical error (see text). The results for the A and SA models are below 0.01, 
independently from X and from (A, A'), and are not displayed. 



It is clear that, if the future experimental results will indicate the LOW solution 
as the preferred one, then anarchical or semi-anarchical schemes will be completely 
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inadequate to describe the data. It would be natural in that case to adopt a model 
where the large gap between the atmospheric and the solar oscillation frequencies 
is built in as a result of a symmetry. What is striking about the results displayed 
in fig. is the ability of the IH models, both in the SS and in the NOSS versions, 
to reproduce the data. An appropriate choice for A (A') is completely successful in 
reproducing r. U e s is numerically close to r and easily respects the present bound. 
Moreover, tan 2 6*12 is very close to 1, due to the pseudo-Dirac structure of the 12 
sector. Only the tan 2 823 distribution shows some flatness and contributes to deplete 
the success rate. It is worth stressing that tan 2 12 is so strongly peaked around 1, 
that any significant deviation from tan 2 8\2 = 1 in the data would provide a severe 
difficulty for the IH schemes. The H model has a smaller success rate (still much 
larger than that of the A and SA models), but has smoother distributions for the 
four observables and is less sensitive than IH to variations of the experimental data. 
The error bars in fig. are dominated by the systematic effects, which have been 
estimated by varying the interval X. We considered four possibilities: X = [0.5, 2] 
(default), X = [0.8, 1.2], X = [0.95, 1.05], and X = [0, 1]. To the highest (lowest) rate 
of each model we then add (subtract) linearly the statistical error. The latter is 
usually smaller than the systematic error. Further details are reported in appendix 
A. 



3.2 LA solution 

The success rates for the LA solutions are displayed in fig. El and El separately for 
the NOSS and SS cases. The two sets of models have been individually normalized 
to give a total rate 100. Before normalization the total success rates for NOSS and 
for SS were in the ratio 1.7:1 (see also appendix A). Although the gaps between 
the rates of different models are reduced compared to the LOW case, nevertheless a 
clear pattern emerges from these figures. The present data are most easily described 
by the IH schemes in their NOSS version. Their performances are better by a factor 
of 10-30 with respect to the last classified, the anarchical models. The ability of the 
IH schemes in describing the data can be appreciated from the distributions of the 
four observables, which, for IH( NO ss,n) an d A = A' = 0.25, are displayed in fig. 01 

The observables r and XJ e z are strongly correlated. Actually, as discussed in ref. 
[T|l2Uj and shown in table 2, in inversely hierarchical models £7 e 3 is typically of order 
r. Therefore, once A and A' have been tuned to fit r, this choice automatically 
provides a good fit to U e3 . Moreover, similarly to the case of the LOW solution, 
tan 2 812 is peaked around 1. At present tan 2 # 12 = 1 is excluded for the LA solution, 
but, thanks to the width of the distribution, the experimentally allowed window is 
sufficiently populated. The width of the distribution is almost entirely dominated 
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Figure 2: Relative success rates for the LA solution, without see-saw. The sum of 
the rates has been normalized to 100. The results correspond to the default choice 
X = [0.5, 2], and to the following values of A = A': 0.2, 0.2, 0.25, 0.3 for the models 
Anoss, SAnoss, IH(noss,ii), and IH(noss,i), respectively (in our notation there are 
no H( NO ss,i), H( NO ss,n) models). The error bars represent the linear sum of the 
systematic error due to the choice of X and the statistical error (see text). 




Figure 3: Relative success rates for the LA solution, with see-saw. The sum of 
the rates has been normalized to 100. The results correspond to the default choice 
X = [0.5,2], and to the following values of A = A': 0.2, 0.3, 0.35, 0.5, 0.15, 0.2 for 
the models A ss , SA SS , H( S s,n), H( S s,i), IH(ss,n) and IH (S s,i), respectively. The error 
bars represent the linear sum of the systematic error due to the choice of X and the 
statistical error (see text). 
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Figure 4: Distributions for IH( N oss,n> 2T = [0.5,2], A = A' = 0.25, obtained with 
10000 points V. 



by the effect coming from the diagonalization of the charged lepton sector. Indeed, 
by turning off the small parameters A and A' in the mass matrix for the charged 
leptons, we get a vanishing success rate for IH( NO ss,i), whereas the rate for IH(noss,ii) 
decreases by more than one order of magnitude. It is worth stressing that even a 
moderate further departure of the window away from tan 2 Q\2 = 1 could drastically 
reduce the success rates of the IH schemes. Finally, the tan 2 623 distribution is 
rather flat, with a moderate peak in the currently favoured interval. All the IH 
models, with or without see-saw, have distributions similar to those shown in fig. HJ 
In particular the tan 2 623 distribution of fig. EJis qualitatively common to all U(l) 
models. This reflects the fact that the large angle 623 is induced by the equal charges 
F(5 2 ) = F(5 3 ), a feature shared by all the models we have considered. To appreciate 
the relevant correlations, the distributions for the model IH( NO ss,ii) are also displayed 
in fig. as scatter plots in the planes (U e 3,r), (tan 2 #i 2 ,r) and (tan 2 623, r). 
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IH(NOSS,II) 




Figure 5: Scatter plots in the planes (U e s,r), (tan 2 9\ 2 , r) and (tan 2 #23,?") for 
IH(noss,ii), 1 = [0-5,2], A = A' = 0.25, obtained with 1000 points V. The box 
shows the experimental window for LA. 

Without see-saw mechanism, the next successful model is the semi-anarchical 
model SA, whose distributions are displayed in fig. |U] (in our notation there are 
no H NO ss,i(ii) models). Compared to the IH case, the tan 2 6 l2 distribution has no 
pronounced peak. Possible shifts in the central value of tan 2 9 12 would not drastically 
modify the results for the SA model. The XJ& distribution is peaked around A = 0.2 
with tails that exceed the present experimental bound. The r distribution is centered 
near r = 1. Finally, the anarchical scheme in its NOSS version is particularly 
disfavoured, due to its tendency to predict r close to 1 and also due to U e z, that 
presents a broad distribution with a preferred value of about 0.5. 

The overall picture changes significantly if the LA solution is realized in the 
context of the see-saw mechanism, as illustrated in fig. El The IH models are still 
rather successful. Compared to the NOSS case, the IHss models slightly prefer higher 
values of r and, due to a smaller A = A', they have very narrow tan 2 9 12 distributions. 
This can be seen from the scatter plots of fig. [7| The other distributions are similar 
to those of figs. 0] and We observe that while most of the points in fig. [7| are 
centered around U e3 ~ 0(A 2 ), there is also a small region clustered at U e3 ~ 0.6. 

Equally good or even better results are obtained by the H(ss,n) model, with 
distributions shown in fig. [S] and EI We see that, at variance with the IH models, 
tan 2 9 12 is not spiky, which results in a better stability of the model against variation 
of the experimental results. The preferred value of r is close to the lower end of the 
experimental window. The U e3 distribution is nicely peaked around A 2 . 
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SA NO ss 




Figure 6: Distributions for the semi-anarchical no-see-saw SA NO ss ; ^ — [0-5,2], 
A = A' = 0.2, obtained with 10000 points V . 



r r r 




Figure 7: Scatter plots in the planes (U e3 , r), (tan 2 12 , r) and (tan 2 #23, r) for IH( SS II ), 
X = [0.5,2], A = A' = 0.15, obtained with 1000 points V. The box shows the 
experimental window for LA. 

The SA S s model is significantly outdistanced from H(ss,ii), IH(ss,ii) an d IH(ss,i)- 
It is particularly penalized by the U e3 distribution, centered around A = 0.3. Finally, 
the least favoured models are H(ss,i) and Ass- The model H(ss,i) fails both in U e s 
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H(ss,ii) 




Figure 9: Scatter plots in the planes (U e3 , r), (tan 2 6*12, r) and (tan 2 6*23, r) for H( SS II ), 
X = [0.5,2], A = A' = 0.35, obtained with 1000 points V. The box shows the 
experimental window for LA. 

(see fig. ITU|) . which tends to be too large for the preferred value of A = A' = 0.5 
and in tan 2 612. The Ass model, as its NOSS version, suffers especially from the U e z 
distribution (see fig. fT0|) which is roughly centered at 0.5, with only few percent 
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of the attempts falling within the present experimental bound. A large U e % can 
be regarded as a specific prediction of anarchy and any possible improvement of 
the bound on \U e %\ will wear away the already limited success rate of the model. 
The distributions in tan 2 6i 2 and tan 2 6*23 are equally broad and peaked around 1. 
Compared to the NOSS case, Ags has a better r distribution, well located inside the 
allowed window. 




As another criterion for evaluating the quality of a model, we address the issue 
of the stability of the observables Oi with respect to small fluctuations of the set of 
coefficients V. Notice that the random coefficients to be put in front of the powers 
of A, A' stand for the combined result of a fundamental theory of flavour, present 
at a certain scale A, and of an evolution from A down to mz- It would thus be 
natural if the physical observables Oi were stable under small perturbations of the 
coefficients, AV, around a given successful representative set Vq. 

This is illustrated for the seesaw case in fig. ^2 where we compare the models 
IH( SS II ), H( SS n) and A ss . The blue dot refers to the observables Oi which follow from 
a typical successful configuration for LA, Vo- The 40 points in red, light blue and 
green are the result of adding to Vo random perturbations AV with \AV\ = \V\j 10 
and random phases. The yellow dots correspond to \AV\ = \V\/2. As appears from 
the scatter plots, IH(ss,ii) is very stable: as already argued in the above discussion, 
actually it is even dangerously stable with respect to the prediction for tan 2 #12. If 
models with hierarchy display a sufficient degree of stability, anarchical ones are 
much less stable, in particular with respect to the predictions for U e3 and tan 2 6* 12 . 
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Figure 11: Results of the stability test in the planes (U e3 ,r), (tan 2 # 12 ,r) and 
(tan 2 # 23 , r) for IH (SS)II) (A = A' = 0.15), H (ss ,ii) (A = A' = 0.35) and A ss 
(A = A' = 0.2), J = [0.5, 2] (see text). 

As shown in fig. H21 the non see-saw case enhances these features: IH( NO ss,ii) is ex- 
tremely stable, while Anoss is higly unstable. For the latter, in particular, even with 
the small fluctuation considered in the figures, tan 2 612 spans from .1 to 10. Thus, 
the criterion of stability supports the same ratings of models previously obtained by 
only considering total rates and distributions. 

As a general comment we observe that our results are rather stable with respect to 
the choice of the interval I. With only one exception, namely the crossing between 
H(ss,i) an d A S s, the relative position of the different models, according to their 
ability in describing the data, does not vary when we shrink X, or when we extend 
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Figure 12: Results of the stability test in the planes (U e3 ,r), (tan 2 6*12, r) and 
(tan 2 6 23 , r) for IH (NO ss,ii) (A = A' = 0.25) and A NOS s (A = A' = 0.2), J = [0.5,2] 
(see text). 

it to cover the full circle of radius 1 3 . This stability would be partially upset by 
restricting to the case of real coefficients V. In that case the relative rates are 
comparable to those obtained in the complex case only for sufficiently wide intervals 
J, typically I« [-1/VX, -\/A]U [VA, 1/VA]. If we further squeeze X around ±1 the 
rates of all models tend to zero. In the real case, the distribution of tan 2 823 is very 
sensitive to the width of the extraction interval X. Indeed, in the flavour symmetry 
basis, charged leptons and light neutrinos mass matrices are both diagonalised by a 
23 mixing angle which tends to 7r/4 when X is squeezed to 1. As an effect, in this 
limit the distribution of 623 is almost empty around 7r/4 and presents two peaks at 
and it/2. On the contrary, in the complex case the distribution of 6*23 is quite 
insensitive to the width of X thanks to the smearing effect of the phases. 



3 The results for X = [0,n] are independent on n, since changing n amounts to perform a 
renormalization of all mass matrices by a common overall scale, which is not felt by the observables 
we have used in our analysis. 
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4 Conclusion 



If a large gap between the solar and atmospheric frequencies for neutrino oscilla- 
tions were finally to be established by experiment then this fact would immediately 
suggest that neutrino masses are hierarchical, similar to quark and charged lepton 
masses. In fact very small mass squared differences among nearly degenerate neu- 
trino states are difficult to obtain and make stable under running in the absence 
of an ad hoc symmetry. Moreover, the presence of very small mixing angles would 
also indicate a hierarchical pattern. At present, the SA solution of solar neutrino 
oscillations is disfavoured so that 2 out of 3 mixing angles appear to be large and 
only one appears to be small, although the present limit is not terribly constraining. 
However, in the case of the LA solution for solar neutrino oscillations, the ratio r of 
the solar to atmospheric frequencies is not so small and it has been suggested that 
possibly this solution does not require a symmetry to be generated. In the anarchi- 
cal framework the smallness of U e3 , assumed not too pronounced, and that of r, are 
considered as accidental (the see-saw mechanism helps in this respect because the 
product of 3 matrices sizeably broadens the r distribution). In the present study 
we have examined in quantitative terms the relative merits of anarchy and of dif- 
ferent implementations of hierarchy in reproducing the observed features of the LA 
solution. For our analysis we have adopted the framework of SU(5)xU(1)f which 
is flexible enough, by suitable choices of the flavour charges, to reproduce all inter- 
esting types of hierarchy and also of anarchy. This framework allows a statistical 
comparison of the different schemes under, as far as possible, homogeneous condi- 
tions. The rating of models in terms of their statistical rates of success is clearly a 
questionable procedure. After all Nature does not proceed at random and a partic- 
ular mass pattern that looks odd could arise due to some deep dynamical reason. 
However, the basis for anarchy as a possible description of the LA solution can only 
be formulated in statistical terms. Therefore it is interesting to compare anarchy 
versus hierarchy on the same grounds. We have considered models both with normal 
and with inverse hierarchies, with and without see-saw, with one or two flavons, and 
compared them with the case of anarchy and of semi-anarchy (models where there is 
no structure in 2-3 but only in 2-3 vs 1). The stability of our results has been tested 
by considering different options for the statistical procedure and also by studying 
the effect on each type of model of small parameter changes. 

Our conclusion is that, for the LA solution, the range of r and the small upper 
limit on U e 3 are sufficiently constraining to make anarchy neatly disfavoured with 
respect to models with built-in hierarchy. If only neutrinos are considered, one 
might counterargue that hierarchical models have at least one more parameter than 
anarchy, in our case the parameter A. However, if one looks at quarks and leptons 
together, as in the GUT models that we consider, then the same parameter that 
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plays the role of an order parameter for the CKM matrix, for example, the Cabibbo 
angle, can be successfully used to reproduce also the LA hierarchy. On the one 
hand, it is interesting that the amount of hierarchy needed for the LA solution is 
just a small power of the Cabibbo angle. On the other hand, if all fermion masses 
and mixings are to be reproduced, a similar parameter is also needed in anarchical 
models, where all the quark and lepton mass structure arises from the charges of 
In comparison to anarchy even the limited amount of structure present in 
semi- anarchical models already improves the performance a lot. And the advantage 
is further increased when more structure is added as in inverse hierarchy models, 
or models with normal hierarchy and automatic suppression of the 23 determinant. 
In the see-saw case all these types of hierarchical models have comparable rates 
of success (except for H(ss,i))- in the non-see-saw versions of inverse hierarchy the 
performance is even better. An experimental criterion that could eventually decide 
between normal and inverse hierarchy models is the closeness of the solar angle 9\2 
to its maximal value. If the data moves away from 7r/4 the probability of inverse 
hierarchy will rapidly drop in comparison to hierarchical models. 
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Appendix: raw data 



In this section we list the results of our numerical analysis, for the case of complex 
random coefficients V. In particular we detail the dependence of the success rates 
on the size of the window X that specifies the absolute value of the coefficients V . In 
most cases the 'systematic' error due to the choice of X is larger than the statistical 
error. The latter is given by 5P/P = \/tJW s , where P = 100 • N s /N is the success 
rate and N s is the number of successes in N trials. We chose N = 10000, 30000, 
50000 respectively, for the LOW solution, for LA without see-saw and for LA with 
see-saw. When fitting the LOW solution, the results for the anarchical and semi- 
anarchical models do not depend on A or A'. 



Table 3: Success rates for the LOW solution, with the see-saw mechanism. 





success rate 


model 


A(= A') 


X= [0.5,2] 


X= [0.8,1.2] 


X= [0.95,1.05] 


X=[0,1] 


A SS 




< 0.01 


< 0.01 


< 0.01 


< 0.01 


SA SS 




< 0.01 


< 0.01 


< 0.01 


< 0.01 


H(ss,ii) 


0.1 


1.7 ±0.1 


2.8 ±0.2 


2.8 ±0.2 


1.6 ±0.1 


H (SS,I) 


0.15 


14.4 ±0.4 


21.6 ±0.5 


23.0 ±0.5 


7.6 ±0.3 


IH(ssji) 


0.03 


41.6 ±0.6 


30.3 ±0.6 


29.0 ±0.5 


31.2 ±0.6 


IH(ss,i) 


0.04 


42.9 ±0.7 


32.8 ±0.6 


31.6 ±0.6 


32.0 ±0.6 
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Table 4: Success rates for the LOW solution, without the see-saw mechanism. 





success rate 


Anoss 


- 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


SA NO ss 




< 0.01 


< 0.01 


< 0.01 


< 0.01 


IH( N OSS,II) 


0.05 


43.4 ±0.7 


33.3 ±0.6 


32.1 ±0.6 


42.9 ±0.7 


IH( NO ss,i) 


0.06 


46.1 ±0.7 


36.5 ±0.6 


35.5 ±0.6 


38.8 ±0.6 



Table 5: Success rates for the LA solution, without see-saw mechanism. 





success rate 


Anoss 


0.2 


0.33 ±0.03 


0.23 ±0.03 


0.25 ±0.03 


0.26 ±0.03 


S A NO ss 


0.2 


2.5 ±0.1 


2.9 ±0.1 


2.6 ±0.1 


2.5 ±0.1 


IH(NOSS,II) 


0.25 


11.4 ±0.2 


8.8 ±0.2 


8.2 ±0.2 


10.2 ±0.2 


IH(NOSS,I) 


0.3 


10.0 ±0.2 


7.9 ±0.2 


7.7 ±0.2 


8.4 ±0.2 



23 



Table 6: Success rates for the LA solution, with see-saw mechanism. 





success rate 


model 


A(= A') 


J = [0.5,2] 


J = [0.8,1.2] 


J = [0.95,1.05] 


J = [0, 1] 


Ass 


0.2 


0.69 ±0.04 


0.62 ±0.04 


0.63 ±0.04 


0.63 ±0.04 


SAss 


0.3 


2.30 ±0.07 


2.12 ±0.07 


2.07 ±0.06 


1.99 ±0.06 


H(ss,n) 


0.35 


4.33 ±0.09 


4.35 ±0.09 


4.34 ±0.09 


2.36 ±0.07 


H(ss,i) 


0.5 


0.58 ±0.03 


0.88 ±0.04 


0.97 ±0.04 


0.27 ±0.02 


IH(ss,ii) 


0.15 


3.92 ±0.09 


3.97 ±0.09 


4.06 ±0.09 


2.57 ±0.07 


IH(ssj) 


0.2 


2.58 ±0.07 


2.24 ±0.07 


2.07 ±0.06 


2.26 ±0.07 
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5 Addendum 



Since the completion of our analysis new important experimental results have been 
published, in particular the first data from KamLAND [21] and the new results from 
SNO [22] • The main new information, of special relevance for the present analysis, 
is that the LOW solution is now discarded and the window corresponding to the LA 
solution has considerably narrowed. In particular the value of r has settled around 
r « 1/35 and the solar angle 6 , 12 is now more than 5a away from maximal. For 
example, in ref. [23] from a comprehensive analysis of all the available data the 
following 3cr window was obtained : 

0.018 < r < 0.053 
\Uesl < 0-23 

0.30 < tan 2 # 12 < 0.64 1 ' 

0.45 < tan 2 ^ 2 3 < 2.57 

These values are to be compared with those in eq. (JBJ), although in eq. (jZJ) we have 
a 3a window according to ref. [23] while the indicative ranges in eq. (JBJ) were less 
precisely defined. These new experimental developments have a large enough impact 
on the present approach to make it certainly worthwhile and timely to update the 
results of our analysis taking the new information into account. As discussed in our 
paper, the deviation of the solar mixing angle from the maximal value is expected 
to produce a pronounced decrease of the success rate of the IH models. Also the 
considerable narrowing of the allowed r range (the upper part of the LA solution is 
now disfavoured) should have important consequences on the relative effectiveness 
of A, SA, H and IH models. It is interesting to precisely quantify these issues and 
this addendum is devoted to this task. 

Given the new values in eq. (J7J), the values of A and A' have been reoptimized 
and in some cases their values are somewhat changed with respect to the previous 
numbers, as shown in Tables 7 and 8, which replace Tables 5 and 6. The values of 
the charges for the representation 10 have also been modified for those cases where 
a rather large value of A follows from the optimization (in order to maintain the 
correct amount of quark and lepton mass hierarchies), as shown in Table 9 that 
replaces Table 1. Since the neutrino mixing parameters are completely independent 
on the 10 charges, this change is only important for a better fit to quark and charged 
lepton masses and mixings. 

The new versions of figs. 2 and 3 that describe the success rates separately for 
the NOSS and SS cases are presented in figs. 13 and 14. The main qualitative 
difference between the new and the old rates is that indeed, both in the NOSS and 
the SS large drop in the IH relative success rate is immediately apparent. 

This is due not only to the smaller upper limit on the solar mixing angle 9 12 but 
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also to the smaller value of r (in fact, we recall from Table 2 that r and 13 are 
expected to be of the same order in IH models). From the updated histograms in 
figs. 13 and 14 we see that normal hierarchy models (with two oppositely charged 
flavons Hn) are neatly preferred over anarchy and inverse hierarchy in the context of 
these SU(5)xU(l) models. In particular, in the SS case, the Hn models with normal 
hierarchy and suppressed 23 sub-determinant are clearly preferred. We recall that 
for the chosen charge values the Hn model is of the lopsided type. 

In conclusion, with all the limitations of the present approach, it is interesting 
that the hierarchical models are preferred over anarchy and inverse hierarchy. In 
particular in the case of see-saw dominance the preferred models are those with 
natural suppression of the 23 sub determinant (e.g those of the lopsided type or 
those with dominance of a light right-handed eigenvalue) which indeed provide a 
simple natural solution of all known constraints. 



Table 7: Success rates for the LA solution, without see-saw mechanism. 





success rate 


model 


A(= A') 


T = [0.5,2] 


J = [0.8,1.2] 


J= [0.95,1.05] 


J =[0,1] 


Anoss 


0.2 


0.02 ±0.01 


0.02 ±0.01 


0.02 ±0.01 


0.03 ±0.01 


SA NO ss 


0.2 


0.44 ±0.03 


0.74 ±0.04 


0.68 ±0.04 


0.45 ±0.03 


IH( N oss,ii) 


0.25 


0.108 ±0.015 


0.106 ±0.015 


0.068 ±0.012 


0.190 ±0.019 


IH(NOSS,I) 


0.25 


0.15 ±0.02 








0.29 ±0.02 
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Table 8: Success rates for the LA solution, with see-saw mechanism. 





success rate 


model 


A(= A') 


J = [0.5,2] 


J = [0.8,1.2] 


1= [0.95,1.05] 


J =[0,1] 


A SS 


0.2 


0.15 ±0.02 


0.11 ±0.02 


0.12 ±0.02 


0.17 ±0.02 


SA SS 


0.25 


0.53 ±0.03 


0.38 ±0.03 


0.39 ±0.03 


0.46 ±0.03 


H(ss,ii) 


0.35 


1.041 ±0.023 


1.027 ±0.032 


0.986 ±0.031 


0.532 ±0.023 


H (SS,I) 


0.45 


0.099 ±0.010 


0.061 ±0.008 


0.016 ±0.004 


0.058 ±0.008 


IH(ss,n) 


0.45 


0.033 ±0.006 


0.007 ±0.003 


0.011 ±0.003 


0.092 ±0.006 


IH(ss,i) 


0.25 


0.020 ±0.003 


0.013 ±0.004 


0.006 ±0.003 


0.035 ±0.004 
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Table 9: Models and their flavour charges. 



Model 


*10 


*5 




(H u , H d ) 


Anarchical (A) 


(3,2,0) 


(0,0,0) 


(0,0,0) 


(0,0) 


Semi-Anarchical (SA) 


(2,1,0) 


(1,0,0) 


(2,1,0) 


(0,0) 


Hierarchical (Hj) 


(6,4,0) 


(2,0,0) 


(1,-1,0) 


(0,0) 


Hierarchical (Hjj) 


(5,3,0) 


(2,0,0) 


(1,-1,0) 


(0,0) 


Inversely Hierarchical (IHj) 


(3,2,0) 


(1,-1,-1) 


(-i,+i,o) 


(o,+i) 


Inversely Hierarchical (IHjj) 


(6,4,0) 


(1,-1,-1) 


(-i,+i,o) 


(o,+i) 
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Figure 13: Relative success rates for the LA solution, without see-saw. The sum 
of the rates has been normalized to 100. The results correspond to the default 
choice J = [0.5,2], and to the following values of A = A': 0.2, 0.2, 0.25, 0.25 for 
the models A NO ss> SA NO ss, IH(noss,ii)> and IH( NO ss,i), respectively (in our notation 
there are no H(noss,i), H(noss,ii) models). The error bars represent the linear sum 
of the systematic error due to the choice of X and the statistical error (see text). 



SNOsalt- LA SS 




Figure 14: Relative success rates for the LA solution, with see-saw. The sum of 
the rates has been normalized to 100. The results correspond to the default choice 
J = [0.5, 2], and to the following values of A = A': 0.2, 0.25, 0.35, 0.45, 0.45, 0.25 for 
the models A ss , SA SS , H (SSjII ), H (SSjI) , IH (SSiII) and IH (SSjI) , respectively. The error 
bars represent the linear sum of the systematic error due to the choice of X and the 
statistical error (see text). 
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